Doubly heavy baryons in QCD sum rules 



Jian-Rong Zhang and Ming-Qiu Huang 
Department of Physics, National University of Defense Technology, Hunan 410073, China 

(Dated: October 30, 2008) 

The mass spectra of doubly heavy baryons are systematically calculated in the framework of QCD 
sum rules. With a tentative heavy-diquark-light-quark configuration, the interpolating currents 
representing the doubly heavy baryons are proposed. Contributions of the operators up to dimension 
six are included in operator product expansion. The numerical results are compatible with other 
theoretical predictions, which may support the (QQ) — (q) structure of doubly heavy baryons. 
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I. INTRODUCTION 

> . 
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The investigation for the doubly heavy baryon containing double heavy-quark Q (Q = c, b) and a single 
light-quark q (q — u, d, s) may provide with valuable information for understanding the nonperturbative 
QCD effects. SELEX Collaboration has reported the first observation of a candidate for a double charm 
i—i, baryon, 1], which was confirmed through the measurement of a different weak decay mode later Q. 
Additionally, BABAR Collaboration and Belle Collaboration have reported that they have not observed 
any evidence of doubly charmed baryons in e + e~ annihilations [3]. The feasibility of doubly heavy baryons 
studied at the Large Hadron Collider (with the design luminosity values of C = 10 34 cm~ 2 s _1 and ^/s = 
14 TeV) was presented 0]. To observe the doubly heavy baryons, it is necessary to give reliable theoretical 
predictions of their properties. So far there have been various approaches, by which the doubly heavy 
t-H ■ baryon masses been calculated. Such as quark models 0, @, 0, S, lil E3| > MIT bag model 11 1 , and mass 
formulas However, the nonperturbative QCD which dominates the low energy physics phenomena 

| has not been fully understood in theory yet. 

Another comprehensive and reliable way for evaluating the nonperturbative effects is the QCD sum rule 
In this method, instead of a model-dependent treatment in terms of constituent quarks, hadrons are 
represented by their interpolating quark currents taken at large virtualities. The correlation function of 
these currents is introduced and treated in the framework of the operator product expansion (OPE), where 
the short and long distance quark-gluon interactions are separated. The former are calculated using QCD 
^ . perturbation theory, whereas the latter are parameterized in terms of vacuum condensates. The result of 
the QCD calculation is then matched, via dispersion relation, to calculate observable characteristics of the 
<_i ' hadronic ground state. Furthermore, the interactions of quark-gluon currents with QCD vacuum fields 



critically depend on the quantum numbers of these currents. However, the accuracy of this method is 
limited, on the one hand, by the approximations in the OPE of the correlation functions and, on the other 
hand, by a very complicated and largely unknown structure of the hadronic dispersion integrals. There 



have the non-relativistic QCD sum rule 14| and the full QCD one [l5| been explored to gain the masses 



of doubly heavy baryons. Nevertheless, Ref. [lfj] has only studied some doubly heavy baryons. Motivated 
by evaluating the mass spectra of them, the analysis of the ground-state doubly heavy baryons via QCD 
sum rules would be done in this work. 

The content of this paper is as follows. In Sec [H] QCD sum rules for the doubly heavy baryons are 
derived. Section [IIII contains numerical analysis and some discussions. 
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II. DOUBLY HEAVY BARYON QCD SUM RULES 

In a tentative picture for doubly heavy baryon QQq system, the low mass q orbits the tightly bound 
QQ pair. The (QQ) — (q) structure may be described similar to Qq mesons, where the QQ pair plays the 
same role of the heavy antiquark Q in Qq. The study of such configuration can help one to adopt the 
appropriate interpolating currents. For the ground states, the currents are correlated with the spin-parity 
quantum numbers + and 1 + for the heavy QQ diquark system, along with the light quark q forming the 
state with J p = A and the pair of degenerate states. For the latter case, the QQ diquark has spin 1, 
and the spin of the third quark is either parallel, J p — | + , or antiparallel, J p — A , to the diquark. The 
choice of and T k matrices in baryonic currents may be determined according to the rules in . For 



the baryon with J p = | , the current may be gained using SU(3) symmetry relations [17J- Consequently, 
the following forms of currents are adopted in present work 

js QQ = e a bc{QlCY k Q b )Y k q c , 

j s . QQ = e abc ^[2{q T a CT k Q b )T' k Q c + {Q T a CY k Q b )Y' k q c ], 

jn QQ = £abc(QaCT 'kQb)Y ' k s c , 

j a . = e abc ^=[2{ s T a CY k Q b )Y k Q c + {Q T a CY k Q b )Y' k s c ], 



js QQ , = £abc{QlCY k Q' b )Y k q c , 

3s- QQ , = e abc -^[(qTCY k Q b )Y' k Q' c + {q T a CY k Q' b )Y' k Q c + {Q T a CY k Q b )Y' k q c ], (1) 

jn QQ , = £abc(QlCY k Q' b )Y k s c , 

m« QQI = £afcc ^|[( s rCT fe Q b )r; ! Q^ + ( s rcr fc Q^r;Q c + (g^cr fc g b )r; Sc ], 

3B> QQ> = e abc (QlCY k Q' b )Y' k q C: 
3n' QQ , = e abc {QlCT k Q' b )Y k s c . 

Here the index T means matrix transposition, C is the charge conjugation matrix, a, b, and c are color 
indices, Q and Q' denote heavy quarks, and q is u or d. The choice of Y k and Y k matrices are listed in 
TABLE |H 

' p 
TABLE I: The choice of F k and T k matrices in baryonic currents. The index d in Sd, Ld, and J d d means diquark. 

{QQ} denotes the diquark in the axial vector state and [QQ] denotes the diquark in the scalar state. 
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The QCD sum rules for the doubly heavy baryons are constructed from the two-point correlator 

U(q 2 ) = i J dW**(0|T[j(aOJ(0)]|0>. (2) 
Lorentz covariance implies that the two-point correlation function has the form 

u( q 2 ) = 4u 1 ( q 2 ) + n 2 ( q 2 ). (3) 

For each invariant function IT and n 2 , a sum rule can be obtained, which is given below. Phenomenolog- 
ically, the correlator can be expressed as a dispersion integral over a physical spectral function 

t 2^ l2 i + M H 1 f°° , ImnP hen (s) 

n (q ) = >?h-T75 ^ + - / ds ^ + subtractions, 4) 

M H -q 2 ir J s-q 2 



where Mh denotes the mass of the doubly heavy baryon. In obtaining the above expression, the Dirac 
and Rarita-Schwinger spinor sum relations, 

J2N(q,s)N(q,s) = tf + M H , (5) 

S 

for spin- 1 baryon, and 

J2 Mq, s)N v (q, s) = {4 + M H )( 9fw \ llllv + q ^\ M ^ - |^), (6) 

for spin-| baryon, have been used. In the OPE side, one works at leading order in a s and considers 
condensates up to dimension six. The s quark is dealt as a light one and the diagrams are considered 
up to order m s . To keep the heavy-quark mass finite, one uses the momentum-space expression for the 
heavy-quark propagator. We follow Refs. [l8[ and calculate the light-quark part of the correlation function 
in the coordinate space, which is then Fourier-transformed to the momentum space in D dimension. The 
resulting light-quark part is combined with the heavy-quark part before it is dimcnsionally regularized at 
D = 4. The correlation function can be written in terms of a dispersion relation as 

f°° P (^) 

Hi(Q 2 ) = / ds— - — -, i = 1,2, m,Q = m,Q> or mg ^ m,Q>, (7) 

J(m Q +m Q ,) 2 S -T 

where the spectral density is given by the imaginary part of the correlation function 

Pi (s) = ilmnp PE ( S ). (8) 

After equating the two expressions for II(q 2 ), assuming quark-hadron duality, and making a Borel trans- 
form, the sum rules can be written as 

X 2 He -M%/M> = / ds Pl { S )e- s l M \ (9) 

J (m Q +m Q ,) 2 

\ 2 H M H e- M */ M2 = f S ° dsp 2 (s)e^ M2 . (10) 

To eliminate the baryon coupling constant A# and extract the Mh, we take the derivative of Eq. §§§ with 
respect to —jji, divide the result by Eq. ([9]) itself, and similarly deal with Eq. ([Toll to yield 

M H = f° ds Pl (s)se- s/M2 1 f° ds Pl (s)e~ s/M2 , (11) 

J(m Q +m Q ,) 2 J(m Q +m Q ,) 2 



Ml = \ dsp 2 { S )se- 8 l M I / ds P2 (s)e- s l M . (12 

J(m Q +m Q ,) 2 J (m Q +m Q ,) 2 

The spectral densities should be distinguished for two kinds of doubly heavy baryons, namely, containin 
the same or differently heavy quarks. Firstly, with 
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III. NUMERICAL ANALYSIS AND DISCUSSIONS 



In this part, the sum rule (fTTj) will be numerically analyzed for better convergence in contrast with 
the sum rule (fT2"|) for this work. The input values are taken as m c = 1.25 GeV,mi, = 4.20 GeV, m s — 
0.13 GeV, (qq) = -(0.23) 3 GeV 3 , (ss) = 0.8 (qq),(gqa ■ Gq) = ml (qq),m 2 = 0.8 GeV 2 , and (g 2 G 2 ) = 
0.88 GeV 4 . Complying with the standard procedure of sum rule analysis, the threshold sq and Borel 
parameter M 2 are varied to find the optimal stability window, in which the perturbative contribution 
should be larger than the condensate contributions while the pole contribution larger than continuum 
contribution. Thus, the regions of thresholds are taken as those presented in the figure captions, with 
M 2 = 3.5 ~ 5.0 GeV 2 for S cc , S* c , Q cc , and fl* cc7 M 2 = 7.5 - 9.0 GeV 2 for S cb , S* fc , fl cb , Q* cb , Z' cb , 
and Q cb , and M 2 = 9.5 ~ 11.0 GeV 2 for S^,, Qtb, and Qt b , respectively. The corresponding Borel 
curves are exhibited in Figs. 1-7. In Table HH the numerical results are collected, together with other 
theoretical predictions. It is worth noting that uncertainty in our results are merely owing to the sum 
rule windows (variation of the threshold sq and Borel parameter M 2 ), not involving the ones from the 
variation of quark masses and QCD parameters. Although the values of masses are in agrement with other 
theoretical predictions, some of the absolute differences from them are not small, for instance, the masses 
of S cc , Q, cc , and S* b , whereas, in the tolerable ranges of the sum rule method accuracy (the upper bound 
of relative accuracy may approximate to 30%). Visually, the Borel curves for 2 CC , f2 cc , and S* fe are not 
very flat, but it is difficult to find much better sum rule windows. That's probably because the condensate 
contributions for them, which may play an important role in stabilizing the Borel curves, nearly vanished 
or are small. The stability of these three curves might be improved by including some higher dimension 
condensate contributions. 
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FIG. 1: The dependence on M 2 for the masses of H cc and H* c from sum rule The continuum thresholds are 

both taken as ,/io = 4.0 ~ 4.2 GeV. 
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FIG. 2: The dependence on M 2 for the masses of Q cc and Q* c from sum rule (JTTJ - The continuum thresholds are 
taken as y/sE = 4.0 ~ 4.2 GeV and y/sE = 4.1 ~ 4.3 GeV. 



-10.6GeV 
-10.7GeV 
10.8GeV 



-10.6GeV 
-10.7GeV 
10.8GeV 



9-5 1<to 2 (GeV 2 ?°- 



5 11 11. 



'« 9-5 1^ 2 (GeV 2 )0-5 '1 11-5 



FIG. 3: The dependence on M 2 for the masses of Hfc& and E£ b from sum rule (JTTJ) . The continuum thresholds are 
both taken as ^/sq = 10.6 ~ 10.8 GeV. 
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FIG. 4: The dependence on M 2 for the masses of Qbb and from sum rule t|lip . The continuum thresholds are 
both taken as Js^ = 10.8 ~ 11.0 GeV. 
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FIG. 5: The dependence on M 2 for the masses of H C 6 and E*;, from sum rule ([lip. The continuum thresholds are 
taken as </sE = 7.4 ~ 7.6 GeV and </sE = 7.6 ~ 7.8 GeV. 
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FIG. 6: The dependence on M 2 for the masses of fl cb and £l* b from sum rule (fTT)l . The continuum thresholds are 
taken as y/sE = 7.3 - 7.5 GeV and ^ = 7.6 - 7.8 GeV. 
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FIG. 7: The dependence on M 2 for the masses of S cb and r2 cb from sum rule (JTTJ . The continuum thresholds are 
taken as ^fsT) = 7.3 ~ 7.5 GeV and 0* = 7.4 ~ 7.6 GeV. 

TABLE II: The mass spectra of doubly heavy baryons (mass in unit of GeV). 
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In summary, the QCD sum rules have been employed to compute the masses of doubly heavy baryons, 
including the contributions of the operators up to dimension six in OPE. The final results are in agreement 
with the existing predictions of other approaches, which may support the heavy-diquark-light-quark con- 
figuration of doubly heavy baryons. Anyhow, the results from all the theoretical approaches are looking 
forward to experimental identification. It is expected that the Large Hadron Collider may supply a gap of 
the doubly heavy baryon data in future. 
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